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NARROW BEAM ANTENNAE 



BACKGROUND OF THE INVENTION 
[0001] Field of the Invention -- This invention 

5 relates generally to antennae and in particular to 
narrow band antennae. 

[0002] The approaches described in this section 
could be pursued, but are not necessarily approaches 
that have been previously conceived or pursued. 

10 Therefore, unless otherwise indicated herein, the 
approaches described in this section are not prior art 
to the claims in this application and are not admitted 
to be prior art by inclusion in this section. 
[0003] To define the direction of radiation D with 

15 angular precision AG requires an aperture of, 

D = 1.22X/Ae (1) 
or so say the text books on optics. Therefore to 
define the direction of a beam in the horizontal plane 
a large area is conventionally needed. 

20 [0004] As illustrated in Fig. 1, devices detect the 

direction of a wave by the oscillations on the surface 
of the detector. A larger detector senses more 
oscillations and is therefore more sensitive to 
direction. 

25 [0005] Basically this is because the wave field has 
to execute a number of oscillations on the 
circumference of this area before we can tell where it 
is coming from: the fewer the oscillations the poorer 
the angular resolution. Mathematically speaking the 

3 0 wave field may be written as, 

Hz=HQGxp{ih'cosQ-i(ot) = HQ ^ J;;, (Ar)/'" exp(//?2e (2) 

m=—oo 

where for illustration a wave polarized with the E 
field in the horizontal plane and the H field parallel 
to the z-axis is assumed. The Bessel function Jj„ is 
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central to the issue of directionality. Roughly 
speaking, 

U/m(^)I~1 kr>m 
« 0 kr <m 

[0006] The number of oscillations of the wave field 
5 around the circumference is restricted by the size of 
kr and hence the limitations on resolution. Fig. 2 
shows a plot of Jjn^{kr) using data taken from M. 

Abramowitz and I. A. Stegun, "Handbook of Mathematical 

Functions" Dover, New York (1972) , which is hereby 
10 incorporated herein by reference in its entirety. 

[0007] Additional references which are useful as 

background to the subject matter contained herein are: 

J. D. Lawson, Journal lEE 95 part III p3 6 3 (1948) ; V. 

G. Veselago, Sov. Phys . USP. 10 509 (1968); J. B. 
15 Pendry, A. J. Holden, W. J. Stewart, I. Youngs, Phys. 

Rev. Lett. 76 4773-6 (1996); J. B. Pendry, A, J. 

Holden, D. J. Robbins, and W. J. Stewart, J". Phys . 

[Condensed Matter] 10 4785-809 (1998); J. B, Pendry, 

A. J. Holden, D. J. Robbins, and W. J. Stewart, IEEE 
2 0 transactions on microwave theory and techniques 47, 

2075-84 (1999); D. R. Smith, W. J. Padilla, D. C. 

Vier, S. C. Nemat -Nasser , S. Schultz. Phys . Rev. Lett. 

84, 4184-4187 (2000); J. B. Pendry, Phys. Rev. Lett. 

85 3966 (2000) ; A, J. Ward, and J. B. Pendry, Journal 
25 of Modern Optics, 43 773-93 (1996); and J. B. Pendry 

and S. A. Ramakrishna, J". Phys. [Condensed Matter] 15 

6345-64 (2003) , all of which are hereby incorporated 

herein by reference in their entirety. 

[0008] As shown in Fig. 2, the Bessel function 
30 7^=9 (/t) controls the amplitude of the 9*^*^ order 
oscillations on the surface of the detector. Evidently 
if the detector is small so that kr«9, then there is 
a weak contribution and the angular sensitivity of the 
device is reduced. 
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[0009] However equation (3) is only approximately 
true. In principle we could take a circle with a 
small circumference and with a highly sensitive piece 
of apparatus measure the amplitude of the high 
5 frequency angular components. This is a severe 
challenge because, although these amplitudes are 
always finite, their magnitude diminishes very rapidly 
as the radius shrinks. To illustrate the point we 
give an approximate expression, 

10 J^{kr)^ ^^ ; kr«m (4) 

[0010] For example for a structure of diameter 
r = X/tt , 

-^".(l)*-^ 1«« (5) 

hence J9 (l) = 2.76x10"^ and this is the magnitude of the 
15 signal we would need to detect for an angular 
resolution of only 3607/w = 40''. The sensitivity 

required increases dramatically as the radius shrinks 
relative to the wavelength. 

[0011] This trade off between sensitivity and 
20 angular resolution for small apertures must always be 
born in mind whatever other means are devised for 
obtaining high resolution. One way or another a very 
sensitive amplifier is required. Conversely, if a 
highly directional signal is to be radiated from a 
25 compact structure, very high intensity fields must be 
injected at some points on the structure. 
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SUMMARY OF THE INVENTION 
[0012] The needs identified in the foregoing, and 
other needs and objects that will become apparent from 
the following description, are achieved in the present 
5 invention, which comprises, in one aspect, an antenna 
comprising a first region having a first refractive 
index, and a second region having a negative 
refractive index, wherein the second region 
substantially surrounds the first region, such that 
10 radiation outside the second region is reproduced in 
the first region, 

[0013] In another aspect of the present invention, 

a method of producing an antenna it taught comprising 
providing a first region having a first refractive 

15 index, and providing a second region having a negative 
refractive index, wherein the second region 
substantially surrounds the first region, such that 
radiation outside the second region is reproduced in 
the first region. 

20 [0014] In the following we present a system for 

creating high angular sensitivity in a compact 
structure. The demands of sensitivity places severe 
demands on the properties and manufacture of the 
components and some of the issues this may involve are 

25 discussed. The central element is the negatively 
refracting materials that have recently appeared onto 
the electromagnetic scene. These have introduced new 
possibilities for control of electromagnetic fields 
and particularly for manipulation of the near fields, 

3 0 which are important ingredients of compact directional 
aerials . 

[0015] To sense the direction of radiation 
precisely with a small diameter aperture the rapidly 
oscillating components that give the directional 
35 information need to be amplified. 
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DESCRIPTION OF THE DRAWINGS 
[0016] The invention will now be described further, 
by way of example only, with reference to the 
accompanying drawings, in which: 
5 [0017] Fig. 1 illustrates how devices detect the 
direction of waves by the oscillations on the surface 
of a detector; 

[0018] Fig. 2 illustrates the Bessel function with 

m=9 ; 

10 [0019] Fig. 3 illustrates a first embodiment of an 
elect romagne tic ant enna ; 

[0020] Fig. 4 illustrates refraction in a negative 
refractive index medium; 

[0021] Fig. 5 illustrates the influence of a 

15 negative refractive index medium; 

[0022] Fig. 6 illustrates an example of a negative 

refractive index material, this material comprising a 
split ring structure; 

[0023] Fig. 7 illustrates a split ring structure 

20 and its permeability; 

[0024] Fig. 8 illustrates the interaction of an 
object on a negative refractive index material; 
[0025] Fig. 9 shows a Cartesian and an cylindrical 
co-ordinate system; 
25 [002 6] Fig. 10 illustrates a wave vector along a 
cylindrical wave guide; 

[0027] Fig. 11 illustrates the objective of a 

narrow beam antenna; 

[0028] Fig. 12 illustrates how co-ordinates of a 

30 cylindrical coordinate system are mapped to planes; 

[0029] Fig. 13 illustrates the variation of 

62 with ^ ; 

[0030] Fig. 14 illustrates the variation of 

i^) with r = yjx^ ^-y^ ; 

35 [0031] Fig. 15 illustrates optical behavior of an 
antenna as illustrated in Fig. 3; 
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[0032] Fig. 16 illustrates the affect on a 
electromagnetic ray of an antenna as illustrated in 
Fig. 3; 

[0033] Fig. 17 is a ray diagram; 
5 [0034] Fig. 18 is a magnetic field diagram of a 
perfect system; 

[0035] Fig. 19 illustrates the magnetic field for 
increasing levels of loss 5; and 

[0036] Fig. 20 illustrates the amplitude d^ of the 

10 mth order of the wave field inside the smallest 
cylinder of an antenna as shown in Fig. 3. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0037] A narrow beam antenna is described. In the 
following description, for the purposes of 
explanation, numerous specific details are set forth 
5 in order to provide a thorough understanding of the 
present invention. It will be apparent, however, to 
one skilled in the art that the present invention may 
be practiced without these specific details. In other 
instances, well-known structures and devices are shown 
10 in block diagram form in order to avoid unnecessarily 
obscuring the present invention. 

[0038] In the following sections, the use of 

negatively refracting materials to compress an 
incoming wave into a smaller volume is discussed. 
15 Fig. 3 illustrates the plan view of an antenna. The 
antenna 20 comprises a first cylinder 3 of radius 
and a second cylinder 2 of radius A2 surrounding the 
first cylinder 3. The outer cylinder 2 comprises a 
negatively refracting material, contained within the 

2 0 cylinder 2 of radius r2 such that, as far as observers 

external to r2 are concerned, it is completely 
invisible. In other words it does not scatter 

incident radiation. The region inside the smallest 
cylinder 3, radius , is filled with a material whose 
25 refractive index is, 

n = ri/r^ (6) 
and within this inner cylinder 3 an observer will see 
a compressed version of the incident wave. The 
compression factor is simply the refractive index, n. 

3 0 This structure maps the contents of a larger cylinder 

1 of radius ri into the smaller cylinder 3, where 

n=nr2=ri/r^ (7) 
[003 9] Fig. 3 shows an embodiment of a narrow beam 

antenna 20. A suitably designed negative material 
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(gray shading) placed in the cylindrical annulus 
between ^2 and will compress the wave field 

originally within the cylinder ri to fit inside the 
smallest cylinder radius . 
5 [0040] We now have a short wavelength version of 
the incident wave travelling in the same direction and 
hence we can employ our detector of choices, e.g. a 
horn antenna or dipole array, to detect the radiation, 
but with the length scale reduced by a factor of n and 

10 therefore with enhanced directionality. It will be 
obvious from reciprocity that an aerial radiating 
short wavelengths inside the small cylinder will 
result in a highly directional beam emerging into 
vacuum. By making a structure of radius we have 

15 gained an effective aperture of radius a-j . The region 
between and r2 is empty space and therefore the 

effective "gain" in aperture is a factor of n^r2/r^ . 
Negative Refraction 

[0041] This simple conclusion results from some 
20 complex mathematics. Negative materials will be 
described in the next section, then the design 
procedure will be explained as well as the materials 
needed to complete the design. 

[0042] The refractive index defines the 
25 relationship between wave vector, k, and frequency, 
(0 , when an electromagnetic wave propagates through a 
material : 

k = rKQ (8) 

where , 

30 n = yfsii (9) 

and e is the electrical permittivity, jx the magnetic 
permeability. We know that if either one of these 
quantities is negative then a wave propagating in such 
a material would result in an imaginary value of n and 
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hence of k . This happens at optical frequencies in 
metals and the imaginary wave vector means that light 
does not penetrate far into a metal and is almost 
completely reflected. 
5 [0043] Some years ago Veselago (V. G. Veselago, 
Sov. Phys. USP. 10 509 (1968)) pointed out that some 
very strange things occur when both s and ix take 

negative values: k is once again real but with a 
strange twist. He argued than whereas we usually 

10 choose n to be positive, in this new situation we are 
forced to choose the negative sign for the square root 
in (9) . Although there has been some heated debate 
about the sign of n the conclusion, now backed by 
several experiments, is that choice of the negative 

15 sign gives the correct results for refraction in 
negative media. 

[0044] Amongst the strange effects noted by 

Veselago was the curious refraction of radiation at a 
surface. The negative refractive index implies that 

2 0 radiation refracts to the "wrong" side of the normal 

giving rise to the chevron style diffraction shown in 
Fig. 4. As shown in Fig. 4, radiation refracts into a 
negative refractive index medium. On the left we see 
the ray diagrams and on the right the direction of the 

25 wave vectors. Note that the wave vector is oppositely 
directed to the group velocity which defines the 
direction of the rays. It follows from these laws of 
refraction that a focussing effect can be achieved by 
a slab of negative material. 

30 [0045] Fig. 5 shows the laws of refraction applied 
to rays emanating from a point source 6 near a 
negative slab 8. Two foci 10, 12 are achieved: one 
(10) inside and one (12) outside the medium 8. In the 
case of n=~l the focussing is free of aberration, but 

3 5 otherwise not so. 
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[0046] As shown in Fig. 5, a negative refractive 
index medium 8 bends light to a negative angle 
relative to the surface normal. Light formerly 
diverging from a point source is set in reverse and 
5 converges back to a point. Released from the medium 
the light reaches a focus for a second time. 
[0047] Materials with 8<0 are relatively easy to 
find. At optical frequencies metals have this 

property, and at lower frequencies a lattice of thin 
10 metallic wires has very similar properties with a 
plasma like form to the dielectric function, 

8 = l-G)^/a)(©+/y) (10) 

[0048] The paper by J. B. Pendry, A. J. Holden, DJ 
Robbins, and W. J. Stewart, IEEE transactions on 

15 microwave theory and techniques 47, 2075-84 (1999) 
shows how to make an artificial material with a 
magnetic response that is effectively negative. This 
is the "split ring" structure shown in Fig. 6. The 
left-hand part of Fig. 6 is a plan view of a split 

2 0 ring 60. The middle part of Fig. 6 shows a sequence 
of split rings shown in their stacking sequence 
separated by a distance i . Each split ring comprises 
two thin sheets of metal. The right-hand view is a 
plan view of a split ring structure 62 in a square 

25 array. Typical dimensions for the split rings are as 
follows : 

Inner radius r = 2 . 0mm 
Width c of each ring = 1 . 0mm 
Spacing d between ring edges = 0.1mm 

30 Lattice constant a,, £ = 10.00mm 

[0049] Fig. 7a shows a split ring structure etched 
into copper circuit board to give negative ]i and Fig. 
7b shows schematic values of the permeability for a 
lossy structure. 



MW/1034342 



-10- 



[0050] Fig. 7a shows one of the early realizations 
of the split rings structures, and typical values 
obtained for the permeability. The figure illustrates 
two important points: firstly negative materials are 
5 strongly dispersive with frequency, often taking a 
strongly resonant form; and secondly loss is often a 
feature in these systems and great care must be taken 
to minimize it. 

[0051] The paper D. R. Smith, W. J. Padilla, D. C. 

10 Vier, S. C. Nemat -Nasser , S. Schultz. Phys . Rev. Lett. 

84 4184-4187 (2000) describes a realization of a 
material with both e<0, |j.<0 and hence with negative 
refractive index. Their pioneering work has been 
confirmed by numerous subsequent studies both 

15 experimental and theoretical. 

[0052] However the key concept for the purposes of 

the present work was introduced in J. B. Pendry, Phys. 
Rev. Lett. 85 3966 (2000) where it was pointed out 
that the focussing action noted by Veselago and 

20 illustrated in Fig. 5 was far more general than had 
been realized. Not only does the slab of n = -l 
material bring the "rays" to a focus, it also acts on 
the near field components of the object forcing them 
to contribute to the image. 

25 [0053] Conventionally the near field dies away 
rapidly with distance and so fails to contribute to a 
conventional image. In contrast the slab of negative 
material actually amplifies the near field and so 
gives the correct contribution to the image of all 

3 0 components, near and far field. This means that the 
image is in principle perfect, through to achieve 
perfection the material must be completely free of any 
loss. However even with a lossy sample it is possible 
to get sub-wavelength resolution and to do better than 

35 a conventional lens. 
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[0054] As shown in Fig. 8, the strongly decaying 
near field of the object excites a surface resonance 
in the negative material and this resonant 
amplification brings the wavefield back to the correct 
5 amplitude. Hence in principle we are able to 
construct a "perfect lens." 

[0055] This ability to manipulate the near field 
with the same precision as the far field is the key to 
designing a compact highly directional aerial. As 

10 mentioned earlier, the high-angular- resolution 
components of the wave field are much reduced in 
amplitude inside a small volume (see Fig. 2) . 
Essentially they become part of the near field and to 
extract the angular information they must be 

15 amplified. This is similar to the problem solved by 
the "perfect lens." To return to Fig. 3, where we 
show the incident wave field compressed into a small 
volume, this compression is obtained by amplifying the 
near field components through resonant excitation of 

20 surface modes in the negative material of cylinder 2. 

There will be very large field intensities present 
within the negative materials and this will be one of 
the issues with which our design must cope . 
Methodology used to Design the Device 

25 [0056] We begin with the philosophy that 

diffraction of waves is easy to understand when the 
geometry is simple. For example refraction of 

radiation at the interface between two dielectrics can 
easily be solved and results in Snell's law for the 

30 change in angle at the interface. Refraction at a 
curved surface is more difficult to calculate 
especially if the radius of curvature is comparable to 
the wavelength. Now we make an analogy: we can 
physically take a slab of dielectric and shape it into 

3 5 a curved surface. For example we could make a 
cylinder from the dielectric. Mathematically we could 
do the same trick by bending the coordinate system 
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changing it from a rectilinear Cartesian system to a 
cylindrical one . 

[0057] Some time ago Ward and Pendry (A. J. Ward, 
and J. B. Pendry, Journal of Modern Optics, 43 773-93 
5 (1996)) showed how to rewrite Maxwell's equations in a 
new coordinate system. They concluded that the 
equations had the same form as the original Maxwell's 
equations, but the bending changed the values of s, pi 
that appeared in the equations. The precise values of 

10 the fields also changed, but what stayed the same were 
the trajectories of rays and the places where they 
came to a focus. We are going to use this idea to 
shape a slab of negatively refracting material into a 
cylindrical lens with the properties outlined. First 

15 we need some mathematics which we take from the 
original Ward and Pendry paper. 

[0058] Fig. 9 shows, on the left, a Cartesian 

coordinate mesh and on the right the mesh for a 
cylindrical coordinate system, 

20 [0059] If Maxwell's equations are rewritten in a 

new coordinate system they take exactly the same form 
as in the old system provided that we re-normalize 8 
and |a according to a simple rule. This affords a huge 
saving in effort because it reduces what appears to be 

25 a new problem to an old one: that of solving for the 
wave field on a uniform mesh in a non uniform medium. 
[0060] Consider a general coordinate transformation 
from a Cartesian, x,y,z , frame to a new set of axes, 

^\ix,y.z), q2{x,y.z). gsi^.y.z). (ii) 

30 [0061] Three unit vectors, Ui,U2,U3, point along 

each of the axes ^i,q^2>^3 respectively. Next we 
introduce the length of a line element, 

2 9 7 (12) 

= 011^^1 ^Qiidqi +Q3zdqi +2(2i2^^i^2 "^^Q^dq^dq^ +2023^^2^% 
where. 
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10 



15 



20 



^ dx dx dy dy 



dz dx 



(13) 



dqi dqj dqi dqj dq^ dqj 

[0062] If the line element is directed along one of 

the three axes qi, 

dSi=Qidqi (14) 

where, 

Qf=Qu (15) 

[0063] Using these quantities we can rewrite 
Maxwell's equations in terms of the coordinates 
9h92y93 ' (see A. J, Ward, and J. B. Pendry, Journal of 
Modern Optics, 43 773-93 (1996) for details). 



(16a) 



(v^xH)'=+8oXe'^'^ (16b) 

which as promised are identical in form to the 
familiar equations written in a Cartesian system of 
coordinates. The re-normalized quantities are. 



b'J = B g'J |ui • (u2 X U3 )| QiQ2Qs (QiQj ) 

A'^" =fx^'^|ui(u2xu3)|0,fi2e3(aey) 



(17a) 
(17b) 



where , 



25 



Uj-U] Ui-U2 Ul'Us 

U2-U1 U2-U2 U2-U3 (18) 
U3U1 U3U2 U3U3 

[0064] The new fields are given by, 

Ej^QjEj^ Hj=QjHj (19) 

[0065] These results greatly simplify if the new 
system of axes is orthogonal, i.e. if ui,U2,U3 are 
orthogonal vectors, since in that case g is a unit 
matrix 
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g = 



1 0 0 
0 1 0 
0 0 1 



(20) 



and. 



Ui-(U2XU3) = 1 (21) 

[0 066] Ward and Pendry tested these transformations 
5 in the case of a cylindrical wave guide with perfect 
metal boundaries and found that the distorted mesh 
shows the same excellent convergence as the original 
uniform version. 

[0067] Fig. 10 shows a three dimensional 

10 implementation of an antenna. In the embodiment shown 
the antenna comprises a first cylinder 3 of a first, 
generally positive, refractive index, and a second 
outer cylinder 2 of a negative refractive index. The 
x-y plane 22 is perpendicular to the longitudinal axis 

15 z of the cylinders. The length of the cylinders 
defines the angular resolution in a plane including 
the longitudinal axis z of the cylinders. 
[0068] In an antenna as illustrated in Figs. 3 and 

10, .an external wave incident on the cylindrical 

2 0 system is compressed into an inner region 3 where it 
has a shorter wavelength and hence can be detected by 
a small device with enhanced angular precision. Fig. 
11 shows the appearance of the system to external and 
to internal observers: to viewers outside the 

25 cylinder 1, the cylinder 1 of radius appears to be 
filled with a radially magnified version of the 
contents of the smaller cylinder 3, radius . Hence 
the system is invisible and appears to be transparent 
to incident radiation. 

30 [0069] Viewed from inside the smaller cylinder 3, 

the world beyond /-3 seems to be filled with a 
miniature version of the region external to ry . To 
observers internal to the circle 3, there is also no 
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evidence of the boundary at r=r;^ which is perfectly 
transparent to outgoing radiation. The material 

between and acts as a wavelength 

expander/contractor. The radial magnification factor 
5 is ri/r3 . Provided that an observer stays outside or 
inside the relevant boundaries, he has no means of 
detecting that the boundaries are there. It should be 
stressed once more that these statements apply only at 
one highly specific frequency at which the device is 

10 designed to operate. 

[0070] Therefore the system is designed which, from 
the point of view of radiation, eliminates the space 
ri>r>r2. The methodology is as follows: to identify 
another system which is known to eliminate a region of 

15 space, possibly of a different shape, and then to 
apply a transformation of coordinates that reshapes 
the known region to the desired region. We shall 
exploit a result recently derived in J. B. Pendry and 
S. A. Ramakrishna, J. Phys . [Condensed Matter] 15 

20 6345-64 (2003) that regions of space filled with 
negative material can optically compensate for 
positive regions: the two regions effectively cancel 
one another from an optical standpoint. For example 
as shown in Fig. 12, if we choose, 

25 e20) = -Bl(2l2-l), A2(l) = -Al(2l2-1), 13<1 <l2 (22) 

[0071] Then not only is the region between I2 and 
I3 perfectly transparent to radiation, but also 
radiation crosses this region with no change of phase 
or amplitude. Optically speaking this space does not 

3 0 exist. The antenna described exploits this result. 

[0072] As shown in Fig. 10, in the new coordinate 
system cylinders are mapped into planes of constant 1 . 
[0073] Since we have shown how to relate coordinate 
transformations to changes in 8,n, we can use the 

35 transformation to deduce what values of 8,|li to choose 
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in the region r'^<r<r^. We divide this region into two 
further domains, 

ei(^) = l^l(/^), r2 <r<ri. 

We then need to determine the values for 81,^.1,82,^2 
5 that the system behaves as required. 

[0074] Having defined the cylinders in the x,y,z 
coordinate we now define a new set of coordinates, 
1,(|),Z so that, in the new coordinate system, the 
boundaries appear not as cylinders but as planes: 

10 x = rQe^/^^ cos^, :)/ = ro^^/^o sin(|), z = Z (24) 

[0075] Next we calculate the quantities needed to 

map from the old to the new coordinate system. From 

equation (13) we have, 

G<t> = ro ^e^^l^ 0 sin^ (j) + e^^l^ 0 cos^ <(> = r^e^l^ 0 
15 [0076] The new s.A follow from (17) , 

and the new fields follow from (19) , 

E\ = Q\ [+E^ cos^+ Ey sin(t)], = \-E^ sm^+Ey coscj)], Ez = QzE^ 
^ (27) 
jfli =Ql[+H^cos^+HySin^'], =g^[-^^sin(|)+/f^cos<l)], Hz=QzH^ 



20 



hence , 

El =ro/loc'/'o[+£^cos<t) + £:^sin<t)], =rQe^^^<'[-E^sm^ + EyCos^'j, Ez = E, 
^1 ='"o/lo^'^'°[+^;ccos(t) + //_^sm(|)], =/-oe'/'o[-//^sm(j) + //_j,cos<|)], Hz=H^ 

[0077] We now come to our result: if in Fig. 12 we 
can define ei,fii and 62, A2 so that they are inverse 
mirror images of one another about the line 



(28) 
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(29) 



(30) 



1 2 = 1 0 l^Crz/'b) then, from an electromagnetic point of 
view, the two regions annihilate one another and, as 
far as the system outside the lines ^3 and (. , is 
concerned, the central region does not exist. We make 
5 our symmetrical choice as follows (it is not unique) . 
[0078] First we set, 

lo = l 

then for the outer region, 
ei = +1, = +1, 

10 for the middle region. 

El =-1, =-1, 

ez =-'-oVl2-21, ln(/-3/ro)<l <ln(r2/ro) 
and for the inner region, 
8, = +1, 8^ = +1, 

8z=+/•oV•^-^•3■^^^ l<ln(/3Ao) 
[0079] We also everywhere set, 
15 ^Al=el, tA<j,=e^, Hz=ez (32) 

[0080] In Fig. 13 we give a schematic plot of 8^ . 
Note the antisymmetry about the point 1 2 = 1 0 ''^('Iz/'b) # 

8(l2+5) = -8(l2-5), 

n(l2+5) = -^(l2-6) ^^^^ 

[0081] This enables us to invoke the theorem that 

2 0 antisymmetrical regions optically annihilate one 
another . 

[0082] Fig. 11 shows the variation of withl , 

Note that within the shaded region 8^ is anti- 
symmetrical about 1 2 = 1 qI^C^'iAg) as required for 
25 focussing- Optically speaking we can remove the shaded 
region and close the gap. 



(31) 
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[0083] Next we use equations (26) to transform to 

the xyz coordinate frame to give for the region 

outside the cylinder radius ^2 , 

^^ = +1, £y=+l, e^=+l, r>r2 

^x = -l. ^y=-h ^z^-r^lr\ r^<r<r2 (34) 

and with identical values for, 



J^x=e^> S^y^^y. V-z^^zl (35) 



[0084] Equations (34) and (35) define the structure 

we wish to create. 

[0085] . Fig. 14 shows a schematic plot of 8^. Note 
10 that 8^ takes the free space value of 8^=1 forr>r2, 
and a constant value for /'<A3- Fig. 12 shows the 

variation of ^zij"^ with radius r = ^/x^Tj^ . Note that 8^ 

is constant outside the "active region" /3<r<r2 where 
the focussing takes place. 

15 [0086] Now let us take the further step of asking 
what happens when we eliminate the gray region from 
Fig. 13. Fig. 15 shows the optical behavior of the 
system taking account of cancellation between the 
antisymmetric regions. As illustrated in Fig. 15a (on 

2 0 the left) , when viewed from a point external to the 
outer cylinder 1, radius ri, the whole of space 
appears as though it were vacuum. As illustrated in 
Fig. 15b (on the right), when viewed from a point 
inside the inner cylinder 3 radius rj the whole space 

25 appears as though ^^V'-r^jr^ . If we view the system 
from outside, then the outer region finishes at, 

ll =loln(ri/ro) (36) 
where the inner region begins. Transforming the 
truncated system back into the xyz coordinate system 

30 leads to Fig. 15a. 
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[0087] Alternat ively if we view the system from 
inside, then the inner region finishes at, 

I3 =loln(A3/ro) (37) 

where the outer region begins. Transforming the 
truncated system back into the xyz coordinate system 
leads to Fig. 15b. Thus, when receiving waves from 
the outside world, the inner region sees the waves as 
obeying Maxwell ' s equations 

' ' ^JiSr^S = ^i'-i'/r^)"} = ±V(^^ A,^)*? (38) 

10 [0088] Since co is fixed and is conserved, the 

in-plane component of wave vector, , is increased by 

a factor / r^; -rx/r^ relative to the incident wave. 

[0089] The maximum aperture of a receiver inside 
the inner cylinder is, 

15 £> = 2r3. (39) 

[0090] However, because the wavelength in the plane 

is reduced by a factor a^/aj , the resolving power is 

equivalent to an aperture in free space of 

^e# = 2r3xri/r3 =2ri (40) 

20 [0091] It is important to note that only the 

region, 

r<r2=yjr^ (41) 
need be filled with material and the region, 

r>r2 (42) 
25 is free space. Therefore the material system is more 
compact than an equivalent free space system by a 
factor of C, where 

C = r2/ri =^3/^2 (43) 
Neqrative Refraction and Curved Surfaces 
30 [0092] In this section we shall explore how the 

structure designed interacts with electromagnetic 
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waves and in particular how the fields are configured 
within the regions where s,n are both negative. We 
begin with a simple ray tracing exercise and ask what 
is the trajectory of a ray which starts at infinity 
5 and interacts with our structure, as shown in Fig. 16. 

[0093] Fig, 13 shows a ray 100 incident on a 
cylinder with impact parameter R . A point along the 
trajectory is defined by the angle (|) , 

[0094] The trajectory has the following form in 
10 free space, 

r = i?/sin<|), r>r2 (44) 

where R is the impact parameter. We also know that 
the trajectory is compressed in the inner region, but 
still executes a straight line because the refractive 
15 index is constant, 

r = Ar/?/sin<|) = -^i?/sin<|), r<ry (45) 

[0095] In the transformed coordinate frame, 

r = roexp(l) (46) 

and therefore the trajectory becomes in the 
2 0 transformed frame, 

/-oexp(l) = i?/sin<j), 1 >l2 (47) 

[0096] Fig. 14 shows a series of ray trajectories 
for the system we have designed. Trajectories 70 
starting at infinity which hit the negatively 

25 refracting cylinder 2 are concentrated within the 
inner cylinder 3. Note the negative angles of 
refraction. In addition there is a set of closed 
trajectories 72 starting within the inner cylinder 3 
which never escape from the system. 

30 [0097] Invoking the symmetry of the trajectory we 

have , 

roexp(2l2 -l) = ^/sin(t), I2 <1 <l3 (48) 
or in the original frame. 
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•j- = R/sin^, r2<r<r2 (49) 

[0098] This is the missing portion of the 

trajectory in the region where the refractive index is 
a function of radius. Remembering that, 

5 ^ = (50) 

we can also write, 

^ = /?/sm(|), r2<r<r^ (51) 
r 

which shows explicitly that (j) has the same value at 
as at . 

10 [0099] Fig. 17 shows a set of trajectories for a 

range of impact parameters. The prediction we made 
that the wave field inside the inner cylinder 3 is a 
compressed version of the contents of the outer 
cylinder 1 were it filled with free space is shown. 

15 The figure has several interesting features that will 
inform our more complete calculation below. 
[0100] First note that only trajectories that 

impact upon the cylinder 2, radius /2 , are refracted 
to the inner cylinder 3 . They are negatively 

20 refracted at both interfaces as required, arrive as 
predicted but fill only the middle half of the inner 
cylinder 3. Only those trajectories which strike the 
middle cylinder 2 are captured, the rest appear to 
escape. This appears to contradict the prediction 

25 that all of the light impacting on the outer cylinder, 
radius ry , will pass through the inner cylinder 3 . 
[0101] The dilemma is resolved by the fact that we 
are dealing with waves and not with rays. The missing 
trajectories are indicated by the numeral 72. We see 

3 0 from the figure that they correspond to closed 
trajectories that never escape from the system. 
However this is true only in the ray approximation. 
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The closed trajectories will in fact correspond to 
cavity resonances and when Maxwell's equations are 
solved will couple to the external rays. Energy will 
leak from the external rays to the internal resonances 
5 and in time the missing trajectories will be populated 
with energy. This process leads to the enhanced 
resolution. 

[0102] Obviously the classical trajectories alone 
contain only the directional information available 

10 from an aperture the size of the middle cylinder 2 
whereas the antenna resolution is of the same order as 
would be obtained in a conventional antenna with an 
aperture equal to that of the outer cylinder radius 
JTj. The effective aperture rj defines the angular 

15 resolution in a plane perpendicular to the 
longitudinal axis of the cylinders. The length of the 
cylinders defines the angular resolution in a plane 
including the longitudinal axis of the cylinders. 
Generally the length of the cylinders is great 

20 compared with the wavelength of the antenna. 

[0103] To calculate how this system responds to an 
incident plane wave of the form, 

-Hz0^^p{^^^os(p-i(ot), r>r2 (52) 
we assume for simplicity that the wave vector is 
25 perpendicular to the axes of the cylinders, and that 
the magnetic field is parallel to the cylinder axis. 
Within the inner cylinder 3 the wave field again has a 
plane wave format but with a shorter wavelength, 

Hzo exp (/^[ri//3 ]r cos cp - /©/), r <r^ ( 53 ) 

30 [0104] Next we must find the fields in the region 

r^<r<r2 which we do by recognizing that 8,n in this 
region have been constructed to be antisymmetric to 
the outer region r2<r<r^. Therefore knowing the 
fields in r2<r<rx we can calculate the fields in 
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r3<r<r2. Substituting from (24) for the new 

coordinates, 1(|)Z : 

H2 = exp^/^Toe^ coscp~/a)/j, 1 > 1 2 ( 54 ) 

Invoking the antisymmetry principle between the two 
5 regions gives, 

= ^20 exp [ikro^e^"^ ^"^ cos cp - /(o/^, 1 3 < 1 < 1 2 (55) 

and substituting back into the original system gives, 

=^zoe^p(''^['2^A]cos(p-/a)/j, A3 <r<r2 (56) 

[0105] Equations (52), (53), (56) specify the 

10 fields everywhere. It is easy to check that these 
fields satisfy the boundary conditions at the two 
interfaces, and are solutions of Maxwell's equations 
in each of their domains. 

[0106] Fig. 15 shows the magnetic field of a plane 

15 wave incident on the system represented as an 
amplitude map. The contours map the phase fronts and 
the wave velocity heads normal to the fronts. When 
compared to Fig. 17 it is noted that the two sets of 
lines are roughly orthogonal as they should be. Note 

20 also that a proper solution of Maxwell's equations has 
now filled in the missing intensity inside the inner 
cylinder 3 where the amplitude is uniform. 
[0107] There are several points of interest in Fig. 

18. Note that the fields are uniform everywhere, with 

25 excursions only between ±1. Contrast Fig. 18 (the 
full solutions of Maxwell's equations) with the ray 
diagram in Fig. 17. The ray picture predicts a region 
of "closed orbits" 72 inaccessible to the outside 
world. However there is no evidence of these empty 

30 regions in Fig. 18: the wave nature of light has 
ensured that these regions are equally populated, at 
least in an ideal system. In this system the ray 
vectors representing the group velocity and Poynting 
vector are orthogonal to the phase fronts, as can be 
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verified by comparing the two figures and remembering 
that in a negative medium the phase and group 
velocities are oppositely directed. 

[0108] The solution to Maxwell's equations in (52), 
5 (53) , and {56) , assumes that the values prescribed for 
8,jJ. can be realized exactly. In the next section we 
examine how the situation degrades when there is less 
than perfect realization of 8,|i and the effects of 
losses which degrade the resonances which populate the 
10 inner cylinder 3 with flux. To calculate the less 
than perfect fields we need to make a more formal 
decomposition of the fields in terms of cylindrical 
harmonics : 

[0109] We assume that region I, r>r2, comprises 
15 free space so that, 

Hi = H'^"" exp(/Vcoscp-/(0/) + iy|''''^^ exp(/A:'i r cos(p*-/co/) 

'"=-^r ml (5 '7) 

= ^zO Z ym{k\r)^<^mH)n\hr)\i^ r>r2 

where , 

^l=<DCo^ (58) 

and Jy^.H^ and cylindrical Bessel functions of order 
20 m. In the absence of any material inside the cylinder, 
there would of course be no scattered wave and a^=0. 
There would also be no scattered wave if the inner 
region 2, 3 was filled precisely as prescribed above, 
but since we are now considering an imperfect 
25 situation, we need to allow for scattering. 

[0110] In region II, r^<r<r2, we no longer have an 

ideal material and therefore writing the solutions of 
Maxwell's equations is more involved. We shall assume 
that the imperfections arise from the materials 
3 0 employed being lossy, and for simplicity we shall 
further assume that the loss takes a particular form: 
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\i{r) = {\^ib)\ii^^^l{r), <r <r2 
e (/-) = (1 + 75) sy^^^/ (r ) , ^3 < r < /2 

[0111] With this assumption we can exploit 
coordinate transformation to write down the solutions 
in this region: 

zO Z [*m-/;;.(^2['-2V'-J) + C;.i';«(^2['-2V'-])|'"exp(,m^ 

where (^2'')>^m (^2^) solve Maxwell's equations in a 
uniform medium with gain, so that, 

k2 =(l-/5)oco^ (61) 

which transforms to absorption when inverted in the 
10 complementary medium. 

[0112] Region III, A3>r, we assume that we can make 

this essentially free of loss since it is a "normal" 
non negative medium. 

/W=+oo 

Hz =^zO 2 ^/w^OT(V)''"exp(/>w(p-/(o/), r3>/- (62) 
m=~-co 

15 where, 

*3 =®Co^nA3 =f^in/^3 (63) 
These assumptions, made to represent imperfections in 
the simplest realistic manner, will not affect the 
qualitative nature of our conclusions . 
20 [0113] It remains to calculate the coefficients, 
which we do in the appendix with the following 
results : 

First we define A^,B^,C^ , 

Am= ^ (64) 

(^2['-2V-3])^:« (V3)-^^^. (V3)^^ 
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o _ Jm{hri)H^m^\hr2)-H^^\k,r2)j^ {k,r2) 

^m- r — (65) 

Jm ( V2)/^£^' {kxr2 ) + -^^i!^ (*1'2 )^ (^2/-2 ) 

YmihriW^ (V2) + ^^i;^(V2)>'«(V2) 

which we then use to calculate the coefficients 
required : 

h -__^L_ r - ^^^^ tc:n\ 
^ _ ^m-^m ( ^2 ) + ^m^";;! (^2^2 ) ~ -^/n ( ) 



(^.[-■jA3])^c„K„(^,[.jA3l) 

"m 77 X (69; 

We now have the solutions we require. 
Discussion of the Materials Required 

10 [0114] Let us recap on our objectives. The aim is 

to achieve angular resolution which beats the 
diffraction limit. The material part of the system is 
contained within r<r2 and the resolution we could 
obtain from a conventional system of these dimensions 

15 is of the order, 

Ae2=V(2^2) (70) 
The fields have been decomposed into cylindrical waves 
of order m: see equation (57) which we reproduce here. 



^i=ffzQ Z + (V) /'"exp(//w(p r>r2 (57) 

20 Clearly the higher values of m contribute greatest to 
the angular resolution and we can give an equivalent 
formula for the resolution, 

^6 = V'"max (71) 
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which from (70) implies that on the boundary /' = r2, 

/?W2 =27cr2A = V2 (72) 

and, 

^im(^1''2)«0> m>kir2 (73) 
5 This is a point made at the beginning in equation (3) • 
[0115] Inside the inner cylinder 3 the fields are 

give by (62) , 

7?I=+O0 

=H^Q dmJm {k^r)f exp{im(p-i(i)t), r^>r ( 62 ) 

[0116] This cylinder is filled with a high 
10 refractive index material so that, 

«3=ri/r3 (74) 

so that, 

h^h^xlr^ (75) 
and in principle this small cylinder 3 has, 

15 ''W3=V3=Vl (76) 

[0117] In other words within this small high 

refractive index cylinder 3 we have the potential for 
angular resolution equal to that provided by a much 
larger cylinder radius in free space. This is only 

20 true if we can ensure that the relevant waves actually 
get through, 

d^^\ m<m^2> (77) 
Our ideal design ensures that d^^^l for all values of 
m but losses will degrade this ideal performance and 

25 hence degrade the resolution. 

[0118] First we discuss the classical ray picture. 
Absorption will reduce the intensity of a ray but 
unless the absorption is very large it will not have a 
catastrophic impact on the ability of a ray to reach 

30 the interior of the inner cylinder 3. Therefore we 
expect the rays mainly to get through. However as we 
have seen in (72) /»niax2 gives only the angular 
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resolution associated with , not the enhanced 
resolution associated with rj , 

[0119] The additional waves with higher values of m 
can only enter the inner cylinder 3 by tunneling 
5 through the resonances corresponding to the closed 
classical ray trajectories. Unfortunately resonances 
are very susceptible to absorption and we expect that 
for even modest values of loss 5 the resonances may be 
sufficiently damped to cut the angular resolution back 
10 to the smaller value. Very low loss materials will be 
needed to realize the full performance of the system. 

[0120] With that introduction let us see what the 
calculations give. We choose a system with dimensions 
appropriate to the GHz regime, 

ri = 0.04m, 
^2 = 0.02m, 
r3=0.01m, 
/ = (o/(27c) = 30GHz 
Thus our ideal resolution would be, 

A01 = X/(2ri ) = 0,01/(2 x 0.04) = 0. 125 radians = 7. 16 degrees (79) 

whereas simply filling the space inside r2 with 
conventional technology would give only half that 
20 resolution, 

A02 = V(2'2) = 0 0l/(2x0.02) = 0.25 radians = 14.32 degrees (80) 
We could be more ambitious choosing a larger ratio 
'iAs / but this places greater demands on the material. 
[0121] We made three sets of calculations 
25 corresponding to, 

5 = 0.1, 

0.01, (81) 
0.001, 

Roughly speaking the larger values of loss are easily 
attained with current technology, whereas 5 = 0.01 will 
be a challenge. 
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[0122] We assume a plane wave incident along the x- 
axis with magnetic field polarized along the z-axis. 
In Fig. 19 we show the real part of H^. Fig. 16 shows 
three separate calculations of for a system 

5 containing a negatively refracting material in the 
middle cylinder 2, and a high refractive index 
material inside the inner cylinder 3 . On the left is 
an overview showing all three cylinders. To the right 
an expanded scale showing just the inner cylinder 3. 
10 The top pair is calculated for low losses, 5 = 0.001, 
the next pair for 5 = 0.01, and the bottom pair for 
5 = 0.1. Note how as the loss is increased the fields 
are confined to the region occupied by the rays in 
Fig. 17. 

15 [0123] Comparing Figs. 18 and 19 the calculations 

show that the very lowest losses of 6 = 0.001 give 
results comparable to the exact calculation. These 
low losses are probably not yet attainable with 
current technology. Higher losses, 5 = 0,01, show the 

2 0 field losing its strength at the top and bottom of the 

green cylinder where the ray trajectories would be 
closed, see Fig. 17. As predicted these regions 
require resonant interactions before the field can 
penetrate which are vulnerable to losses. Finally for 
25 5 = 0.1 we see considerable departures from the ideal 
fields. Inside the inner cylinder 3 the fields are 
more or less confined to the open ray trajectories. 
Outside the outer cylinder 2 there is evidence of 
strong forward scattering from the system, and overall 

3 0 the field intensity is much reduced even where the 

rays penetrate. Our qualitative conclusion from Fig. 
19 is that losses of the order of 6 = 0.01 or less are 
desirable . 

[0124] Fig. 20 shows the amplitude of the mth 

3 5 order component of the wave field inside the smallest 
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cylinder 3 (see equation (62) ) for various levels of 
loss as measured by 5 . The more non zero values of 
the better the angular resolution. The vertical 
line to the right shows the number of non zero values 
5 we are aiming for and would give the maximum 
resolution, the vertical line to the left shows what 
conventional technology would give by exploiting the 
space inside r2 . 

[0125] To estimate the angular resolution we can 

10 expect in each case we show in Fig. 20 the cylindrical 
wave amplitudes within the smallest cylinder 3, d^^ , 
see equation (62) . The maximum non zero value 

dictates the angular resolution. The vertical marker 
to the right shows Wmaxl the ideal we are aiming for. 

15 The extremely low loss calculation, 5 = 0.000001, 
comfortably exceeds this target, but with this set up 
the resolution is in any case limited by Wmaxl • The 
marker to the left shows /Wmaxl which is attainable by 
conventional technology. Therefore a cut off in d^ 

2 0 above ^max2 means that an improvement in the 
diffraction limit has been achieved. 

[0126] The high loss material shows no improvement 

over conventional technology. Intermediate losses of 
5 = 0.01 give useful improvement from A62 =14.32'' to 
25 around AG = 11°, low loss of 5 = 0.001 gives A9 = 9°, 
approaching the ideal of A0=7°- 

[0127] Reference has been made above to an antenna 
of a generally cylindrical form. However, the antenna 
may be provided in any suitable form. For instance 
30 the antenna may comprise, as illustrated in Fig. 21, 
an inner sphere 103 of a first refractive index 
substantially enclosed within an outer sphere 102 of a 
negative refractive index. Other geometrical forms 
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are envisaged and the above description relating to a 
cylindrical form is not intended to be limiting on the 
implementation of an antenna as taught . 

[0128] The spherical equivalence of equations (34) 
5 and (35) are as follows: 

^2 

Bj,=By =8^=+-^, 0<r<r^ 
^3 
^2 

8^=8^=8^=-^, r^<r<r2 (82) 
r 

8^ = 8^ = 8^ = -hi, r2 < r < 00 

\^x=^x> V^y=By. V^z^^z (83) 

[0129] Thus there is provided an antenna which 

beats the diffraction limit for the angular 

10 sensitivity of an aerial of a given aperture. The 
device exploits negatively refracting materials to 
enhance the angular resolution. In an ideal situation 
infinite improvement is in principle possible. A test 
calculation which aims to improve resolution by a 

15 factor of two, and assumes that ideal materials were 
available, shows that the design is a successful one. 
Consideration of losses which are likely to occur in 
real materials reduces this improvement but for a 
material having a level of losses which we believe may 

20 be attainable in the near future resolution could be 
substantially enhanced relative to the conventional 
limit. The materials challenge is severe, but when 
low loss negative materials become readily 
commercially available, the diffraction limit on 

25 angular resolution may be beaten effectively. 

[013 0] Although the foregoing description of the 
present invention has been shown and described with 
reference to particular embodiments and applications 
thereof, it has been presented for purposes of 

30 illustration and description and is not intended to be 
exhaustive or to limit the invention to the particular 
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embodiments and applications disclosed. It will be 
apparent to those having ordinary skill in the art 
that a number of changes, modifications, variations, 
or alterations to the invention as described herein 
5 may be made, none of which depart from the spirit or 
scope of the present invention. The particular 
embodiments and applications were chosen and described 
to provide the best illustration of the principles of 
the invention and its practical application to thereby 

10 enable one of ordinary skill in the art to utilize the 
invention in various embodiments and with various 
modifications as are suited to the particular use 
contemplated. All such changes, modifications, 

variations, and alterations should therefore be seen 

15 as being within the scope of the present invention as 
determined by the appended claims when interpreted in 
accordance with the breadth to which they are fairly, 
legally, and equitably entitled. 
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